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ABSTRACT 

Combining Ha and IRAC images of the nearby spiral galaxy NGC 628, we find that between 30- 
43% of its 8/im dust emission is not related to recent star formation. Contributions from dust heated 
by young stars are separated by identifying H II regions in the Ha map and using these areas as a 
mask to determine the 8/xm dust emission that must be due to heating by older stars. Corrections 
are made for sub-detection-threshold Hn regions, photons escaping from Hn regions and for young 
stars not directly associated to Hn regions (i.e. 10-100 Myr old stars). A simple model confirms this 
amount of 8/im emission can be expected given dust and PAH absorption cross-sections, a realistic 
star-formation history, and the observed optical extinction values. A Fourier power spectrum analysis 
indicates that the 8/xm dust emission is more diffuse than the Ha emission (and similar to observed 
Hi), supporting our analysis that much of the 8/im-emitting dust is heated by older stars. The 8/mi 
dust-to-Ha emission ratio declines with galactocentric radius both within and outside of Hn regions, 
probably due to a radial increase in disk transparency. In the course of this work, we have also found 
that intrinsic diffuse Ha fractions may be lower than previously thought in galaxies, if the differential 
extinction between Hn regions and diffuse regions is taken into account. 



1. INTRODUCTION 

The 8/im dust emission is an attractive star forma- 
tion rate measure for galaxies, as it provides excellent 
resolution in nearby low-redshift galaxies. It is also ac- 
cessible at higher redshift as the rest-frame 8/im shifts 
into other accessible infrared bandpasses. Indeed, the 
mid-infrared dust emission around 8/im broadly corre- 
lates with star formation rate indicators such as Ha 
(IRoussel et al.ll200lt IPeeters et al.ll200i IWu et al.1120051: 
IKennicutt et al.l 12009ft. Paa (ICalzetti et al.l 120051 12007ft . 
far inf rared (IPeeters et al J 12004 IDale et all 12005ft and 
radio (iVogler et all l2005TlWu et all 12005ft. However, 
ICalzetti et al.l (|2005l ) draw attention to the non-linearity 
of the 8/im dust emission with respect to the number of 
ionizing photons from Hn regions, raising some concerns 
about its use a s a sta r formation indicator. Furthermore, 
ICalzetti et al.l (|2005ft also note that the 8/im emission is 
more diffuse than either recombination line emission from 
ionized gas (Ha, Paschen-a) or 24/im emission tracing 
the hot dust localized to star-forming regions. The goal 
of this paper is to quantify the galaxy-wide contribution 
to 8/im dust emission from non-star forming sources as a 
guide for using such emission as a star formation tracer 
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in both local and high-redshift galaxies. 

8/im dust emission is made up of two main compo- 
nents, aromatic band emission and continuum emission 
from warm dust grains. In general, the aromatic bands 
lie between 3-17/iin, with some of the strongest bands at 
3.3, 6.2, 7.7, 8.6, 11.2, and 12.7 /mi. Each of these bands 
has been identified with vibrational mode s of C-C and 
C-H bonds of atoms in aromatic rings (see lTielensi r2008. 
for a recent review). Polycyclic aromatic hydrocarbons 
(PAHs) are now the widely-accepted carr iers of the mid- 
IR ba nd f eatures (originally sugges ted bv lLeger fc Pugetl 
119841 and lAllamandola et all 11985ft . although other hy- 
drocarbon compounds may contribute as well. For sim- 
plicity, we will attribute the emission features to PAHs 
for the remainder of the paper. Because PAH emis- 
sion is very bright in photo-dissociated regions (PDRs) 
surrounding sites of active star formation, it was ini- 
tially surmised that PAH emission might be a very good 
star formation rate indicator, at least on global scales 
(|Roussel et al.|[200lft . 

The second component to 8/im dust emission is contin- 
uum emission. Large dust grains may emit at 8/im if they 
are heated by a very intense radiation field (significant 
when the radiation field density is > 10 4 times local). 
But much more important for the 8/im emission is the 
contribution from very small (PAH-sized; « 6A) grains 
that are stochastically heated to high effective temper- 
atures. The importance of single-photon heating was 
made clear by Infra Red Astronomical Satellite (IRAS) 
observations of diffuse regions which show much more 
emission at 12 and 25/im than expected for classically- 
sized grains, given the low radiation field inte nsity 
(jBoulanger et al.l Il985h iBoulanger fc Peraultl 11988ft . To 
fit these and other observations, both PAHs and very 
small grains are now included in dus t grain models (e.g. 
iDesert et ai1ll990HLi fe Drainell2001ft . 

The PAH band emission is now known to have a 
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cirrus component. Band emission associated with the 
diffuse interstellar medium (ISM) in the Galaxy was 
first i nferred from photomet ric measureme nts from the 
IRAS (jBoulanger et al.l[l985h and AROME (IGiard et alJ 
1994) and is easily visible in the Spitzer IRAC 8/xm im - 
ages from the GLIMPSE survey (jBeniamin et al.ll2003l ). 
Spectroscopic confirmation that the photometric excess 
(over the expected dust emission) is indeed due to 
the aromatic emission ban ds came from the I nfra Red 
Telescope in Space (IRTS; iTanaka et al.l 119961) and the 
Infrared Space O bservatory (ISO; iMattila et al.1 119961 : 
ISakon et al.ll2004[). Using mid- IR spectra of interarm re- 
gions. iVogler et al.1 (|2005l) and ISakon et al.1 (|200n ) iden- 
tify the aromatic bands in the diffuse ISM in the spiral 
galaxies M83 and NGC 6946, respectively, using ISO and 
AKARI. Interestingly, these spectroscopic data show dis- 
tinct variations of the PAH bands (for example, relatively 
stronger in the 11.3/im band) in the diffuse r egions, when 
comp ared to activ ely st ar-forming regions ()Vogler et al.l 
120051: ISakon et al.ll2004 120071) . This band variation sig- 
nals a difference in the PAH molecule population and/or 
ionization. 

Using Spitzer IRS spectrosc opy on the SINGS s ample 
of local star-forming galaxies, Smi th et al.1 ()2007l ) sepa- 
rate the PAH bands from the rest of the 8/im dust emis- 
sion, establishing that typically 80% of the 8/im dust 
emission is from the PAH bands (their Figure 12, bottom 
panel) , with presumably the remainder due to continuum 
emission from warm, stochastically-heated grains. How- 
ever, scatter from galaxy to galaxy clearly exists and low- 
metallicity galaxies in particular tend to have low ratios 
of PAH band emission to warm dust continuu m emission 
(|Engelbracht et al .1120051: 1 Cannon et al.ll2006D . 

In this paper, we set out to measure the fraction of 8/im 
dust emission that is not due to recent star formation. 
While the line dividing 'recent' star formation from older 
generations of stars is somewhat arbitrary, here we will 
count populations of stars under 100 Myr. The goal is 
to provide a quantitative warning about using the 8/tm 
dust emission from a galaxy as a direct tracer of star 
formation. We take the nearby face-on galaxy NGC 628 
as a test case and tie this determination to Ha emission, 
which we assume is entirely powered by photoionisation 
from the youngest stars. In Section 2, we present the 
Ha and Spitzer IRAC data used. Section 3 describes the 
methods used to first separate 8/tm dust emission related 
to H II regions and then emission related to stars under 
10 Myr, stars between 10-100 Myr and stars older than 
100 Myr. In Section 4, we present and discuss our results, 
while in Section 5, we test various systematics within our 
method. Section 6 presents our conclusions. 

2. DATA 

2.1. NGC 628 

We chose to work on the unbarred Sc galaxy NGC 628 
becau se it is relatively nearby (7.2 Mpc: iKennicutt et al.1 
1201 lh . large in angular size (10.5 x 9.5 arcmin 2 ), virtu- 
ally face on (i=6°) and observed in all the bands (Ha, 
8/xm, 3.6/tm) we require. NGC 628 is not classified as an 
AGN either by optical em ission lines dHo et al.lll997f) o r 
by mid-IR emission lines (|Goulding k, Alexander! 120091) . 
Its stellar mass is 3.6 x 10 9 M (|Skibba et al.l 120 111) 
and its current star formation rate (SFR) is 0.68 M 



vr~ 1 ()Kennicutt et al.l 1201 lh . Based on spectral energy 
distribution fitting. lAniano et al.l (|2012l ) determine that 
~11.6% of NGC 628's infrared luminosity comes from 
dust heated at high radiation field intensities (U > 10 2 ), 
a typical value among the SINGS sample spirals. 

2.2. Ha data 

We use Ha data obtained at the Kitt Peak National 
Observatory (KPNO) 0.9 m telescope i n Jan uary 1997 
and originally presented in iGreenawaltl (fl998) . Narrow- 
band line and broad-band continuum exposures were 
taken with 16 dithered exposures of 900s and 480s, re- 
spectively. The narrow-band filter used is narrow enough 
(27 A) to exclude contamination from the nearby [N n] 
line. 

A Ha line-only image was obtained after subtracting 
a scaled broad-band image with the scale factor deter- 
mined by comparing the broad-band to narrow-band flux 
ratio of 94 foreground stars. Using the IDL MMM pro- 
cedure on these flux ratio values, we reject outliers and 
determine a mean of 0.806 with a standard deviation of 
0.036 (84 non-rejected stars). The standard error on the 
sample mean is thus only 0.0039 (i.e. only 0.5%); how- 
ever, this assumes that the foreground stars correctly 
approximate the narrow-band to broad-band flux ratio 
of NGC 628's average stellar population. In light of this 
possible systematic effect, we will test values changing 
within the full la range of stellar values in Section 5.3 
(thus ±0.036 or 4.4%). 

Flux calibration was based on standard star observa- 
tions and the conversion to emission measure performed 
using the formula appropriate for Ha and a T — 10, 000K 
plasma: 
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The point spread function (PSF) of the final Ha image 
has a full-width half-maximum (FWHM) of 1.8 arcsec 
and the images have an rms noise level of 3.8 pc cm' 6 
with a pixel size of 0.68 arcsec 2 . 

2.3. IRAC data 

The 3.6 and 8 /tm data on NGC 628 were obtained 
from the SINGS Data Release 5. A tilted plane was fit 
to both images to subtract the background, after using 
an iterativ e algorithm to determ ine the set of background 
pixels, see lAniano et al.l (|2012) for full details. 

Next, we created a 8/tm dust-only map by subtract- 
ing a fraction of the /„(3.6 /tm) image (having first con- 
volved it to the 8 /tm resolution) in order to remove stel- 
lar emission from the /„ (8 (tm ) image. We subtr acted 
0.255 times 7„(3.6 /mi), as in faelou et all (|2004l) . but 
also check a range of believable subtraction fractions in 
Section 5.2. We also note that subtracting the 3.6 /tm 
band will remove slightly more than the stellar emission 
as the 3.3 /xm PAH band is within the 3.6 /tm IRAC 
filter. However, the 3.3 /im feature should be only 0.03- 
0.09 times the strength of the 8/tm feature (jFlagev et al.1 
2006) which should be around 0.8 of the 8/im dust emis- 
sion. So using a factor of 0.255 for the stellar subtraction, 
we will only cause a 0.6-1.8% over-subtraction of the dust 
emission at 8 /tm. 
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The rms noise level in the non-stellar j//„(8/im) im- 
age (with a pixel size of 0.75 arcsec 2 ) is 1.3 X 10" 10 
erg s " 1 cm" 2 sr" 1 and the measured PSF FWHM is 
2.82" (|Aniano et al.|[201lh . 

3. ANALYSIS 

In order to separate the 8/im dust emission linked to 
recent star format ion, we first define Hii re gions using 
the hiiphot code (|Thilker et all 120001 120021 ) on the Ha 
data. H II regions signal the presence of massive, ionizing 
stars, which must have ages under 10 Myr. By masking 
the emission from these regions (see right-hand panel of 
Fig. [J), we calculate the amount of 'non-H II region' emis- 
sion for both the Ha and non-stellar J„(8/im) images. 

This method has two caveats. First, UV photons 
can escape from Hn regions, heating PAHs outside of 
their defined boundaries. Secondly, we miss some star- 
formation by only using Ha emission. Some low mass 
young clusters will be under our detection threshold and 
so included in the non-H II region emission. Similarly, a 
minority of H n regions maybe very heavily attenuated 
and thus not detectable in Ha. Below, we explain how 
we put limits on these contributions, thereby arriving at 
an amount of 8/Ltm dust heated by >10 Myr stars. 

A further consideration is that stars between 10 and 
100 Myr should contribute signi ficantly to the he ating of 
PAHs outside of Hn regions (jPeeters et al.ll2004T) , These 
stars are related to 'recent' star formation, as the term 
is usually used by the extragalactic community, so we 
will separately estimate the contribution of these stars 
to the non-stellar 7„(8/im) image, finally arriving at a 
firm lower limit to how much 8/im dust emission must 
be related to purely old stellar populations. 

3.1. Identification of HII regions 

hiiphot was run on the optical data with default pa- 
rameters, including a signal-to-noise (S/N) cutoff of 5. 
Different emission measure (EM; units of pc cm -6 ) slope 
cutoff values of 1, 1.5, 2, and 4 cm" 6 were specified. 
These different AEM cutoffs control how far H n regions 
grow from initially identified Hn region 'seeds'; smaller 
values of AEM result in larger H n regions. A AEM cut- 
off of 1.5 cm" 6 is the recommendation from lThilker et all 
(2000), a value motivated by wanting to recover as much 
flux spatially linked to an ionizing sources as possible. 
But we also keep maps based on a wider range of AEM 
cutoff in order to investigate how the definition of an H II 
region changes our results (see Section 5.6). 

The noise level and resolution of the data affect how 
hiiphot finds H II regions and we investigate these sys- 
tematics in Section 5.4. Because we wish to define H II re- 
gions on a scale relevant to the 8/xm image with its larger 
(2.8") PSF, the minimum resolution we can achieve with 
this analysis is a physical scale of 98 pc. Thus we con- 
volve the Ha image to reach this resolution and also re- 
sample to a pixel size of 1.32" = 43.6 pc pixel" 1 . In the 
convolved Ha image, the noise is at a level of 3.8 pc cm" 6 
(la). 

The output of hiiphot produces (among other prod- 
ucts) an integer mask identifying H n region pixels (val- 
ues > 1), non-H II region pixels (0), and pixels masked 
from the start (-32768; bright foreground stars). Fig. [T] 
shows the outlines of the determined H n regions (with 
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Figure 1. Hn regions as determined by HIIPHOT in magenta con- 
tours over the Ha image (upper) and stellar-continuum subtracted 
8 fim image (lower). Grey circles indicate regions with foreground 
star contamination, masked out during the analysis. 

AEM=1.5 cm" 6 ) on both the Ha image and over the 8 
/zm dust image. 

3.2. Attenuation correction 

The Ha map must be corrected for Ha photons lost 
to dust absorption and scattering. Integral field spectro- 
scopic data of NGC 628 determines V-band attenuations 
(Ay) using the Balmer decrement method wherever both 
Ha a nd H/3 emission line s are detected (mostly Hll re- 
gions; [SlncheiIeF^l][20Tl|) . No radial variation in the Ay 
values of these regions is found, with an average Ay of 
1.24 (but large scatter at all radii). C omparing Ha and 
24/im emission, Presco tt et al.l (|2007| ) calculate Ay val- 
ues in 500 pc apertures centered on bright 24 /im peaks, 
similarly selecting for Hn regions. Using this method, 
a radial gradient is found, with A# Q declining by -0.53 
dex out to i?25- Figure [2] shows the comparison of the 
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Figure 2. The three different extinction corrections used as a 
function of radi us. The solid bl ack line is the Hn region cor- 
rection from Prcscott ct al. (2007) (updated to lower values from 
ICalzetti et all 11200711 " calibration coefficient and the offset in Ha 
image calibration), while the solid light blue line shows a con- 
stant H II region corr ection factor corresponding to 1.24 in Ay from 
[Sanchez ct al. ( 2 0111 ). The dashed red line gives the linearly de- 
clining extinction correction for the non-H II regions derived from 
a dust surface density map with an added cont ribution from the 
Milky Way extinction from (Schlcgel e t alj|1998ft . 



two H n-region extinctions, together with that for diffuse 
regions (see below) . As both the average gas density and 
the metallicity decline with radius, we select the radially- 
declining version of the H n region extinction correction 
as the better p hysically motivated. H owever, we update 
the values from iPrescott et all ([2007), correcting the ef- 
fect of an incorrectly calibrated Ha image and using the 
newer calibration coefficient of C alzetti et al.1 ([20071 ) to 
derive a central A# Q of 1.03. 

Evaluating the attenuation in the diffuse regions of 
galactic disks is even more difficult. He re, we choose 
to use the dust surface density map from lAniano et al.l 
(2012) in order to estimate the attenuation in the diffuse 
regions. Matching the Hn region attenuation, we assume 
a linear functional form. We hnd a dust surface density 
of 2 x 10 5 M Q kpc~ 2 in a diffuse region toward the center 
of the galaxy and 5 x 10 4 M Q kpc -2 at the edge of the 
detected region in the dust m ap, at approximately 9 kpc. 
Converting to extinction via (|Aniano et~al1l2012f i: 



A v [mag] = 0.67 



1 x 10 5 Mq kpc' 



(2) 



we obtain the Ha attenuation, after additionally convert- 
ing from V band to Ha. However, we assume that the 
emission, on average, originates in the midplane of the 
galaxy, so we take half of these attenuations, as shown 
by the red dashed line in Fig. O Additio nally, we apply 
the fo reground Milky Way extinction from lSchlegel et all 
(p9l . 

3.3. Annular photometry 

We choose to analyze the data in large annular re- 
gions, in order to have enough signal-to-noise to de- 
termine the non-H n region emission level even in the 
outer regions of NGC 628. The annuli all have width 
of 0.1i?25 = 31.4"=1.1 kpc. Using the integer masks 
produced by hiiphot, we compute the Hn region and 
non-H II region flux in each annular region (hereafter de- 
noted i^Ha) and F(8fim), where F(8/im) is defined by 
vFy). Then, to account for diffuse emission contained 
within the Hn regions but not actually related to the 



Hn regions, we take the average non-H n region flux per 
pixel within each annulus and multiply by the number 
of Hn region pixels. This value is subtracted off the 
annular H n region flux and added to the annular non- 
Hn region flux for both Ha and 8/xm data. While the 
/„(8/im) image is already background subtracted as ex- 
plained in Section 2.3, a large (width of 50"=1.8 kpc) 
outer annulus around the whole galaxy is used for back- 
ground subtraction of the Ha image from 100 to 150" 
beyond i?25- 

Using this method and not yet applying an extinction 
correction, we determine a non-H II region F(Ra) frac- 
tion of 37% (commonly called a 'diffuse' Ha fraction). 
We note t hat this value is low er (by 7%) than that re- 
ported in iThilker et al.l (|2002f) , who use identical data 
and the same code. This difference results from the dif- 
ferent accounting of diffuse flux within the Hn regions. 
Here, we have subtracted the average non-H II region sur- 
face brightness from within a given Hn region's radial 
annulus, while IThilker et al.l (|2002[ ) use a more sophis- 
ticated approach that interpolates the background using 
information immediately outside each H II region. As the 
Ha surface brightness is generally higher immediately 
outside of Hn regions, this approach results in higher 
diffuse levels subtracted for each Hn region. Here, we 
would only like to subtract the Ha emission essentially 
from 'in front' of the Hn region, not that linked to it, 
hence taking the annular average approach. 

3.4. F(8fim) contributions from < 10 Myr stars 

The non-H n region 8/im flux cannot be directly 
equated with the 8/im flux produced by stars older than 
10 Myr for two reasons. First, photons escape from H II 
regions and secondly, low-mass or highly-embedded clus- 
ters may not be detectable in the Ha image. 

Unfortunately, integrated UV escape fractions for H II 
regions have not received much attention in the literature 
(line-of-sight attenuations are more the norm) , although 
escape fractions due to ionizing (Lyman continuum) pho- 
tons have. Such studies generally conclude that leaking 
ionizing photons power the diffuse Ha emission, since 
other potential sources (shocks from mechanical feedback 
or in-falling gas, photoionization from post-AGB stars 
and X-ray binarie s) appear unable to provide enough 
ionizing flux (e.g. Fer guson et al.l [19961 ) . Indeed, spec- 
troscopic data confirm that photoionization is do minant, 
altho ugh not the only ionizing source in galaxies ([Martini 
119971 ). How the ionizing radiation permeates so far away 
from Hn regions remains partly mysterious (e.g. iSeonl 
2009) ; but the inhomogencity of the ISM and the low re- 
combination rat e in low density reg i ons probably suffice 
as explanations dMiller fc Coxl[l993l: [Dove fc Shulll ll994l: 
iZurita et all l2002f ). iDong fc Draine) ([20111 ) argue that 
much of the diffuse Ha emission may come from regions 
not in ionization equilibrium, including emission from 
plasma recombining after the photoionization source is 
removed. 

In any case, escaping Lyman continuum photons are 
almost certainly accompanied by longer- wavelength UV 
photons, the photons most effective at exciting PAHs. 
While the interception of ionizing versus UV photons is 
due to different sources (i.e. atomic hydrogen and dust, 
respectively) , if the escape fraction is primarily mediated 
by physical holes in the cocoon of dust and gas around 
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Hn regions, the escape fractions should be similar. Al- 
though the exact ratio between UV and ionizing photon 
escape fractions is not well constrained, we will assume 
the fractions are equivalent. Later, we also discuss the re- 
sults if twice the fraction of UV photons escape compared 
to the ionizing photon fraction (see Section 5.1). Thus, 
the non-Hll region F(Ha) also gives a signal about how 
much UV flux has escaped H II regions and is available to 
excite PAHs and heat small dust grains. To estimate the 
amount of 8/im emission powered by these escaping UV 
photons, we take the diffuse Ha emission and multiply 
by the F (8/im) /F (Ha) ratio found for the Hn regions. 

One caveat to using the non-Hll region F(Ka) to cor- 
rect for escaping UV photons is that some of the non- 
H II region Ha is scattered off of dust instead of signal- 
ing in-situ ionization. Recent studies of high-latitude 
diffuse regions in the Milky Way estimate that approx- 
imately 20% of the Ha emission is originally from disk 
H II r egions and scattered by dust back down to wards the 
disk (fWitt et al.ll20Tol ; iBrandt fc Draind [20121) . A simi- 
lar fraction is likely to hold in the non-Hll regions and 
we apply such a correction here, subtracting 20% from 
the non-Hll region F(Ha) before using it to correct for 
escaping UV photons from Hn regions. 

Secondly, we consider the contribution from very young 
stars (<10 Myr) that lie outside of our defined Hn re- 
gions. Our minimal Ha luminosity region detected is 
« 5 x 10 36 erg s -1 . An individual 07 V star produces 
an Ha luminosity of ~ 10 37 erg s _1 , so we can detect 
regions that only host one or very few ionizing stars. 
However, some low-mass clusters will not have any ion- 
izing stars and thus not be detected using hiiphot. We 
find that H n regions with low Ha luminosity have lower 
F(8[im) I F(Ha) ratios (by about a factor 2) than more 
luminous H n regions. This smaller ratio indicates that 
while these low Ha luminosity clusters have an ioniz- 
ing star or two, their IMF is not fully sampled and they 
have fewer than average bright-UV emitting stars (B- 
type). Assuming other, equally young, clusters have 
higher B-star to ionizing star ratios, we would also like to 
count their contribution to the <10 Myr F(8fim). The 
twice lower F(8fim)/F(Ha) indicates that we are missing 
about half of the 8/im emission that would be measured 
with a fully sampled IMF. The total Ha luminosity of 
the low-luminosity (but detected) clusters, is only 5% of 
the total, so we estimate we need to re-attribute about 
2.5% of the non-Hll region F(8/im) to be associated to 
<10 Myr stars. 

We also likely miss some very highly attenuated H II re- 
gions. However. iPrescott et al l ()2007l ) compared Ha and 
24/im maps for the SINGS galaxies and determined very 
few (< 4%) such regions exist in normal spiral galaxies. 
As a rough correction, we subtract 4% of the Hn region 
-F(8/im) from the non-Hll region F(8^im) and attribute 
it to H II region _F(8/im). Combining the sub- ionizing and 
embedded contributions, the subtraction factor becomes 
6.5%. 

We first perform this correction for sub-ionizing young 
clusters (Equations 3 and 4), then the correction for es- 
caping UV photons (Equation 5 and 6): 



■F(8//m) nonH ii 

= F(8/iTO) nonH ii,i 



0.065 x F(8/im) H n,, 



F(8/iTO) H ii = 1.065 x F(8/im) H n,, 



(3) 



(4) 



F(8/im)>i0Myr 

= F(8/iTO) nonH n - F(Ha) nonHU F -}ftt m ) mi , (5) 

F(Ha)mi 

F(8/im)<ioMyr = F(8/im)tot - F(8^m) >WM yr- (6) 

3.5. F(8[im) contributions from 10 to 100 Myr stars 

Stars remain UV bright and thus effective at heating 
PAHs and small grains up to approximately 100 Myr. 
Such young stars are still tied to recent star formation, 
so we would like to estimate their contribution to the 
8/im dust emission. For both the PAH and dust heat- 
ing, we estimate the contribution based upon their ex- 
pected UV flux (912-3000 A) of 10-100 Myr stars. Using 
iBruzual fc Chariot! (|2003t ) models, with a solar metallic- 
ity and Chabrier IMF, we find that the expected UV flux 
of stars between 10-100 Myr is 0.49 times that of the stars 
younger than 10 Myr. Note that these models are based 
on an exponenti ally decaying SFR with a timescale of 
5 Gyr (eqn. 6 of INoeske et aLl[2007h . 

Assuming that the UV photons of 10-100 Myr stars are 
equally effective at producing 8/im dust emission, then 
0.49 times the F(8/im) from < 10 Myr stars should be 
produced by 10-100 Myr stars. So, we have: 

-F(8/Wl)<i00Myr = 1-49 X F(8/fm)<ioMyr- (7) 

However, this almost certainly leads to an overestimate 
because the 10-100 Myr stars are probably less effective 
at heating PAHs and dust for several reasons. First, they 
are generally further away from dust clouds so the num- 
ber of their photons intercepted by dust is likely lower. 
Secondly, the older stars' UV emission is not as blue and 
thus will have a reduced cross-section. Thirdly, some of 
the 8/im dust emission from H n regions is actually a sec- 
ondary product of the Lyman continuum photons, with 
PAHs excited by Lyman-a, Balmer lines, and 2-photon 
continuum. As we are most concerned with providing 
a lower limit to the amount of F(8/im) from >100 Myr 
stars, we keep the overestimate as a safe approach, thus: 



F(8/i77i) ldst a 



> F(8/im) tota l - 1-49 X F(8/iTO)<i Myi- 

'(8) 



4. RESULTS 



4.1. Integrated values 

Excluding star- forming regions identified by hiiphot, 
we arrive at a total (out to R25) non-HII region dust 
F(8/im) fraction of 71% in NGC 628. Then, correcting 
for other contributions from very young (<10 Myr) stars 
we arrive at a fraction of 62%. Considering the amount 
of UV produced by 10-100 Myr stars (see Section 3.5), 
we estimate a lower limit to the amount of 8/im dust 
emission heated by sources other than young stars as 
43%. These are all much higher values than the non-H 11 
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Figure 3. F(8fJ,xa) fractions and F(8/J,m)/ F(Ha) ratios as a func- 
tion of radius. The top panel shows the non-HII region F(8fim) 
(gray), the >10 Myr F(8/im) (red), and the lower bound for the 
>100 Myr F(8/im) (orange). The attenuation-corrected non-Hn 
region F(Ha) fraction is also shown (blue). The total percentages 
of each quantity are shown on right-hand side of the plot. The 
bottom panel shows the F(8/J,m)/ F(Ha) ratio as it declines with 
radius for the non-Hn (solid line) and Hn (dotted line) regions. 

region F(Ha) fraction, determined to be only 24%, after 
applying the differential extinction correction. 

Thus the non-star forming component of the 8/xm dust 
emission is significant compared to the emission directly 
associated with young stars in NGC 628. Star formation 
rates computed based on the total 8/im dust emission 
thus trace partly the current star formation rate tied to 
the O and B stars, but also take into account emission 
that is heated by the integrated population of stars found 
in the field. This latter contribution may change from 
galaxy to galaxy and thus pollute the estimates of star 
formation rates. Future work on additional galaxies will 
show if this is the case. 

4.2. Radial trends 

Fig.[3]shows the radial trends of the dust F(8/im) frac- 
tions contained outside of Hn regions (grey line), heated 
by stars older than 10 Myr (red dotted line) and heated 
by stars older than 100 Myr (orange dashed line). The 
blue dot-dashed line shows the fraction of Ha emission 
located outside of Hn regions. 

For all of the F(8/im) fractions, we see a dip with a 
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Figure 4. F(8iim) / F(Ha) ratio for individual Hn regions as a 
function of their galactocentric radius. As with the annular aver- 
age, a clear decline is seen with increasing radius. No relation with 
HII region size (shown by symbol color) is seen. 



minimum at 4.4 kpc, the radius at which the bright- 
est Hn regions are found in NGC 628. With the 8/im 
dust emission dominated by these bright regions in this 
annulus, a lower fraction of non-star formation related 
8/im dust emission is expected. After this dip, a grad- 
ual increase with radius is seen. This increase is proba- 
bly due to the relatively greater amount of quiescent gas 
(and dust) at larger radii. Thus as star forming regions 
become more sparse in the outer regions of the galaxy, 
photons from older stars become correspondingly more 
dominant in producing 8/im dust emission. 

We also note the decline of the F(8/mi)/ F(Ha) ra- 
tio for both the Hn regions and the non-Hn regions, 
seen in the bottom panel of Fig. [3] This decline is ex- 
pected because the gas density declines with radius and 
thus the disk transparency (assuming a relatively con- 
stant dust-to-gas ratio) increases. There are simply fewer 
PAH molecules to absorb stellar photons, while there are 
still sufficient hydrogen atoms to absorb hv > 13.6 eV 
photons. Potentially, there is also a contribution from 
the metallicity gradient within the disk as PAH s are ob- 
served to be deficient at low metallicities (e.g. | Maddenl 
l2000HEngelbracht et al.ll2005HGalliano et al.ll2005D This 
observed radial trend highlights the importance of per- 
forming this analysis in radial bins. 

4.3. F(8fim)/F(Ha) ratios in individual Hn regions 

Fig. [4] shows the F (8/im) /F (Ha) ratios of individual 
Hn regions against their galactocentric radius. To make 
this plot, only HII regions with 8/im surface brightnesses 
higher than 10% above the subtracted background were 
used (832 of 1055 regions), avoiding those regions whose 
individual 8/im flux is not well determined. As previ- 
ously seen with the annular mean, the ratio clearly de- 
clines with radius, but there is a very large spread at 
all radii (up to a factor 100). Neither Hn region size 
(shown as symbol color in Fig. |4|) nor Ha luminosity (not 
shown, but similar appearance) appear to correlate with 
the F(8/mi)/F(Ha) ratio. 

4.4. Fourier power spectrum analysis 
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Figure 5. Fourier power spectrum of the Ha (blue) and I u (8fim) 
dust (red) images. At large scales, the 8(im dust emission has 
a clearly steeper slope than the He* emission. The steeper slope 
signals more relative power at large scales, confirming that more 
of the 8^tm dust emission is spread on larger scales than the Ha 
emission. As we are primarily concerned with the slopes, we have 
not applied any particular normalization to the power spectra. 



As an independent test, we perform a Fourier power 
spectrum comparison between the Ha and /„(8/im) im- 
ages. We use the native-resolution and original pixel 
scale images, deprojected to a face-on orientation and ig- 
noring scales below twice the PSF FWHM transformed 
into the deprojected plane. We use IDL to compute 
the 2-dimensional discrete Fourier transform and then 
take the square of this complex array to arrive at the 
power spectrum image. We azimuthally average this 2- 
dimensional power spectrum to obtain the power at each 
linear scale. This 1-dimensional average power spectrum 
is shown in Fig.[5]for both Ha and /„(8//m) dust images. 

The slope of the Ha power spectrum is clearly shal- 
lower than that of the 8/xm dust emission. We choose to 
fit the power spectrum slope from 20-200 arcseconds (w 
0.7-7 kpc) in both images, as both appear approximately 
power-law in this region. The power-law slopes are —0.9 
and —1.5 for Ha and I„(8/im) dust images, shown as 
solid lines in Fig. [5j respectively. The value found for the 
8/xm dust image is comparable to the —1.6 slope found 
for the H I e missio n over similar (0.8-8 kpc) scales by 
iDutta et al.l ([2008). This similarity signals that there 
is indeed 8/im dust emission distributed like the neutral 
(and not necessarily star-forming) gas. 

The —0.9 power-law slope found for the Ha image 
is somewhat similar to two studies that find slopes of 
—0.7 and —0.8 f or the large-scale p ower- la w slope for Ha 
images of M33 lElmegreen et al.l (|2003[) : iCombes et all 
(2012). However, not all galaxies have such shal low Ha 
slopes, in the same paper. lElmegreen et al.l (|2003l ) report 
a slope of —1.5 for NGC 5055 which is closer to the pre- 
diction for a turbulent medium (predicted value of —1.7). 

In NGC 628, the relative steepness of the dust I„(8/im) 
power spectrum compared to that of the Ha image shows 
that the 8/im dust emission is distributed on larger scales 
than the Ha emission. This result fits nicely with the 
larger non-H n region dust F(8[im) fraction than for the 



non-H II region .F'(Ha) fraction we find in the above anal- 
ysis. 

4.5. A simple model 

It has long been suggested that old stars significantly 
contribute to the FIR emission of galaxies, both from 
studies of our own Galax y and e xt ernal galaxies (e.g. 
Lonsdale Persson fc Heloul Il987t iBuat fc Deharvengi 



19881: iTrinchieri et all 119891: iBloemen et al.l 11990: 
Sodroski et all 119971; iDraine et al.l I2007T ). Furthermore , 
Groves et al.l (|2012| ) have recently demonstrated that 



exclusively old stars heat the dust found within the 
bulge of M31. Here, we demonstrate the plausibility 
that stars older than 100 Myr produce a significant 
fraction of PAH emission in NGC 628, while the more 
detailed calculations appear in the Appendix. 

We model the stellar emiss ion of NGC 628 using the 
evolutionary synthesis code of lBruzual fc Char iot (2003) 
and four possible star formation histories (constant, 
linearly decreasing, exponentially decreasing, peaked 
at 2 Gyr). The parameters for the linearly de- 
creasing and 2 Gyr-peaked mod els are based u pon 
iFraternali fc Tomassettil (|2012D and lZou et all (|2011| ) , re- 
spectively, while a decay time of 5 Gyr is used for the 
exponential model. For each stellar age bin, we then 
calculate the absorption of this starlight by applying 
the du st absorption cross-section from the IDraine fc Lil 
(2007) Milky Way dust model with a PAH mass fraction 
of 4.6%. The scaling of the dust optical depth is tied 
to the optical attenuation, as determined in Section 3.2. 
The dust geometry is assumed to be a shell for the very 
youngest stellar populations, fading into a mixed star 
and dust geometry for the older regions, although we 
note that uniformly applying a mixed geometry produced 
nearly identical results. The fraction of dust absorption 
due to PAHs at each wavelength is separated by compar- 
ing the PAH and total dust cross-sections. This allows 
us to capture the effect of increased UV cross-section of 
the PAHs compared to larger dust grains. 

Figure [6] shows the fractional contributions of different 
age bins for two different star formation histories (con- 
stant and exponentially declining) to the total stellar, 
stellar UV, TIR and PAH emission for the model (atten- 
uation values fixed to those at the center of NGC 628). 
These plots show that the fractional PAH luminosity 
from the young stellar population is indeed higher than 
that of the total TIR, due to the stronger UV absorption 
of these small molecules. However, both have signifi- 
cant emission attributed to heating from stars older than 
100 Myr as shown by the fractions shown for the lowest 
bar of each plot of Fig. [51 26% and 35%. While these frac- 
tions of old-star heated PAH emission are smaller than 
what we derive observationally for NGC 628 (43%), shift- 
ing to a reasonable twice-as-high escape fraction for the 
UV as for the ionizing photons lowers the observed frac- 
tion to 30% (see Section 5.1), similar to the results of 
the model. The model and observations appear to be 
consistent within the uncertainties. 

5. TESTS OF SYSTEMATICS 

5.1. UV escape fraction 

Arguments can be made as to why the UV escape frac- 
tion should be similar to that of the ionizing photons, but 
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Figure 6. The fractional contribution of each age bin to the to- 
tal stellar, FUV (912-3000A) stellar, TIR and PAH luminosity. 
The age bins for the young populations (< 100 Myr) are shown 
in shades of blue and each bin represents 20 Myr. The age bins 
for the old populations are shown in shades of red and each bin 
represents 2 Gyr (other than the first, which is 1.9 Gyr). The top 
four bars represent a constant SFH stellar population while the 
bottom four bars represent an exponential stellar population with 
timescale 5 Gyr. 



because different physical mechanisms are in play, the ra- 
tio between the fractions is unlikely to be unity, as we 
have assumed. In fact, the observed diffuse -F(Ha) frac- 
tion is only 28% while the escape fraction calculated for 
the typical Hn-region Ay of 0.7 at 1500 A is 49% (ap- 
plyin g the starburst attenuation law from ICalzetti et al.l 
1994). It is thus likely that the UV escape fraction is 
larger than that for the ionizing photons. Assuming twice 
the rate of escape for UV photons compared to ionizing 
photons changes dust F(8/im) fraction from 62% to 53% 
for >10Myr and from 43% to 30% for >100Myr stars. 
This change obviously lowers the F(8fim) attributed to 
old stars, to a level more in line with our simple model. 
Our observations thus suggest that the old-star powered 
8/zm dust emission fraction lies within these ranges. 

5.2. Stellar subtraction at 8[im 

A fraction (ris^m) of the 3.6 /mi IRAC band intensity, 
I„(3.6 /im), is subtracted from the I v (8 /tm) image in 
order to obtain a stellar-subtracted 8 /jm image. We test 
values of ry^m from 0.22-0.30, centered on the adopted 



value of 0.255 of lHelou etaTI (|2004ft . Non-H n region dust 
F(8[im) fluxes are about 3% greater in the ?78 M m=0.22 
images than the ry8/im=0.30 images, while the Hn region 
dust F(8/xm) are only about 1% greater. As both of these 
quantities simultaneously increase, the non-Hn region 
dust F(8fjxa) fraction only ends up changing by 0.5%, a 
negligible amount. 

The non-Hn and Hn region F(8/xm) dust/Ha ratios 
change by the full 3% and 1%, respectively. But these 
few-percent changes are small compared to radial decline 
seen in both the non-H n and H n region ratios. 

5.3. Ha stellar subtraction 

Similarly, a fraction of the wide-band image (mostly 
stellar continuum) is subtracted from the original 
narrow-band Ha image (containing both line and contin- 
uum emission) in order to obtain an Ha line-only image. 
Here, we test the effect of this stellar subtraction on our 
results. In Section 2.2, we noted that the mean ratio 
between the narrow-band and broad-band fluxes for 84 
foreground stars was 0.806 with a standard deviation of 
0.036. We use the ±1<7 values to determine how different 
our results could be with different ?7hq- 

The subtraction affects the non-H n regions (-10%, 
15%) more than the bright Hn regions (-2.5%, 1%) and 
especially in the galaxy center where the starlight is most 
concentrated. But based upon inspection, the different 
Ha subtractions do not significantly affect the spatial ex- 
tent or number of regions selected by hiiphot. Overall, 
the non-H II region F(Ha) fraction changes from 34% to 
40% for the under- and over-subtracted images (36% for 
best-value) if no extinction correction is applied (to ei- 
ther Hn or diffuse regions), or only from 22% to 27% 
(24% for best-value) if our preferred extinction correc- 
tion is applied. So, as with the 8/im stellar subtraction 
experiment, it appears that the choice of Ha stellar sub- 
traction value does not strongly influence our results. 

5.4. Ha image properties 

Two properties of the Ha image (resolution and noise) 
influence how hiiphot identifies HII regions. The res- 
olution here was set to match that of the 8/im images, 
at 2.82" = 98 pc at the 7.2 Mpc distance of NCC 628. 
Reducing the resolution to 4.23"(148 pc, in both the Ha 
and 8/xm images) increases the diffuse F(Ha) fraction 
from 24% to 27% and thus ends up decreasing the frac- 
tion of F(8fj,m) heated by old stars from 47% to 41%. 
Testing the effect of using a finer scale would be more 
informative, but the fact that there is such a decrease 
as we degrade the resolution signals that even less of the 
Ha emission may truly lie outside of H n regions and thus 
even more of the 8/im emission is in fact heated by old 
stars. Thus our conclusion that this number is a lower 
limit will hold. 

We add Poisson noise in both the narrowband and 
broadband images separately and re-subtract to obtain 
a line-only image. We add more and more noise (mim- 
icking a proportional decrease in exposure time between 
narrow and broad band images) until a rms noise level 
of 20 pc cm -6 is obtained in the line-only image. Us- 
ing such noisier data does not change any of our results 
within a 2% margin. 

Increasing the noise by a significant factor was a less 
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important effect, producing only a 2% change in the final 
estimate of the fraction of F(8/im) powered by old stars. 

5.5. Dust attenuation effects 

Correcting the Ha emission for dust attenuation sig- 
nificantly changes the non-H n region (i.e. 'diffuse') frac- 
tion of F(Ha). This effect can be seen clearly by fol- 
lowing the blue line in the upper panels of Fig. [71 With 
absolutely no extinction correction applied, the non-H n 
F(Ha ) fraction is 3 7 %. W ith the Hn region correction 
from iSanchez et al.l (1201 If) and the non-H II region cor- 
rection from the lAniano et al.l (|2012[ ) dust map, the frac- 
tion is 25%. With o ur adopted ext inction correction o f 
iPrescott etafl (|2007| ) (corrected via lCalzetti et al.ll2007l ) 
for H II regions and again the dust map for the non-H II 
regions, it is only 24%. Thus, considering the differential 
extinction between within and outside of H II regions is 
very important to determining the intrinsic amount of 
diffuse Ha emission. Values may therefore be lower than 
the 50% commonly reported by studies that do not 
take such differential extinction into account. 

As the amount of diffuse Ha emission determines how 
much diffuse 8/im we attribute to emission escaping from 
Hn regions and thus associated to young stars, the ex- 
tinction correction noticeably propagates through to our 
8/im fractions attributed to >10 Myr and >100 Myr pop- 
ulations, as seen in Fig. [7J Accounting for differential 
extinction between the H n regions and diffuse regions 
leaves about 10% more 8/im emission to be attributed 
to old stellar populations. Thus the extinction correc- 
tion is an important systematic, although the difference 
between our two adopted extinction corrections is fairly 
minimal, only at the 2% level. 

5.6. Emission measure cutoff effects 

The emission measure cutoff specifies the gradient in 
EM (measured in units of cm -6 pc) at which the growth 
of Hn regions is terminated by hiiphot. Lower EM cut- 
offs result in larger Hn regions. Figure [3] shows that 
decreasing the EM cutoff lowers both the F(Ha) and 
F(8fim) non-H II fractions, as expected. These changes 
are significant, with the non-H II F(Ha) fraction decreas- 
ing from 30% to 20% (Prescott extinction correction ap- 
plied), the non-H II F(8/mi) fraction from 76% to 66% 
and the lower limit on the F(8fim) from old stars from 
49% to 38%. However, inspection shows that the Hn re- 
gions for AEM of 4 misses much of their contiguous Ha 
emission, and those for AEM of 1 are likely too large. 
Still, we expect the uncertainty on defining the exact 
boundary on Hn regions to contribute an «5% uncer- 
tainty on our derived lower limit to the amount of 8/im 
dust emission that is heated by stars older than 100 Myr. 

6. CONCLUSIONS 

Masking out H n regions with the help of an Ha image 
and the HIIPHOT code, we find that 30-43% of the non- 
stellar F(8fim) in NGC 628 is due to PAH molecules and 
dust heated by stars older than 100 Myr. Our analysis 
shows that 38-47% is heated by stars under 10 Myr, while 
21-23% is heated by stars between 10 and 100 Myr. The 
major parameter contributing to the ranges of these es- 
timates is the unknown escape fraction of UV photons 
compared to ionizing photons. 



We thoroughly investigate the systematic effects of 
background subtraction, stellar emission subtraction (for 
both 8/im and Ha images), extinction correction and 
emission measure cutoff and estimate a combined sys- 
tematic uncertainty of «10% on the lower limit for non- 
star formation related F(8^im). Thus this work rigor- 
ously establishes that a significant fraction of F(8fim) in 
NGC 628 is not due to heating from the youngest popu- 
lation of stars. The quantity of this non-star formation 
related emission may be different from galaxy to galaxy, 
complicating measures of the star formation rate deter- 
mined using 8/im emission at both high and low redshift. 

We also find that the F(8fim)/F(H.a) ratio declines 
with galactocentric radius in Hn regions by a factor 5.4 
and outside of Hn regions by a factor 4.0. This decline 
is likely caused by the easier destruction or more difficult 
formation of PAHs in low metallicity regimes. However, 
PAH abundance may not be the full answer. Changes 
in PAH population, in terms of size distribution or ion- 
ization state, could also cause a reduction in 8/im band 
emission. 

In the process of this work, we find that the fraction 
of non-H II region (commonly called diffuse) Ha emission 
is lower than usually reported, taking into account the 
higher preferential extinction of Hn regions. With our 
adopted extinction model, we find an intrinsic non-H II 
region F(Ha) fraction of only 24%. This may be com- 
pared to 37% when differential extinction between HII 
regions and diffuse emission is not taken into account. 

Independent of the Hn region based analysis, we com- 
pare the power spectra of the 8/im dust and Ha images 
and find that the 8/im dust image has more power at 
lower spatial frequencies (larger scales) compared to the 
Ha image. This distribution confirms the more diffuse 
nature of the 8/im dust emission. The power-law slope 
of the 8/im dust power spectrum matches well with that 
of the Hi emission in NGC 628, suggesting that the 8/im 
emission is more strongly linked to the Hi gas than to 
star formation. 

REFERENCES 

Allamandola, L. J., Tielens, A. G. G. M., & Barker, J. R. 1985, 
ApJ, 290, L25 

Aniano, G., Draine, B. T., Gordon, K. D., & Sandstrom, K. 2011, 

PASP, 123, 1218 
Aniano, G. et al. 2012, ApJ, 756, 138 
Benjamin, R. A. et al. 2003, PASP, 115, 953 
Bloemen, J. B. G. M., Deul, E. R., & Thaddeus, P. 1990, A&A, 

233, 437 

Boulanger, F., Baud, B., & van Albada, G. D. 1985, A&A, 144, 
L9 

Boulanger, F. & Perault, M. 1988, ApJ, 330, 964 
Brandt, T. D. & Draine, B. T. 2012, ApJ, 744, 129 
Bruzual, G. & Chariot, S. 2003, MNRAS, 344, 1000 
Buat, V. & Deharveng, J. M. 1988, A&A, 195, 60 
Calzetti, D., Kinney, A. L., & Storchi-Bergmann, T. 1994, ApJ, 
429, 582 

Calzetti, D. et al. 2005, ApJ, 633, 871 
— . 2007, ApJ, 666, 870 
Cannon, J. M. et al. 2006, ApJ, 647, 293 
Churchwell, E. et al. 2006, ApJ, 649, 759 
Combes, F. et al. 2012, A&A, 539, A67 

da Cunha, E., Chariot, S., & Elbaz, D. 2008, MNRAS, 388, 1595 
Dale, D. A. et al. 2005, ApJ, 633, 857 

Desert, F.-X., Boulanger, F., & Puget, J. L. 1990, A&A, 237, 215 
Dong, R. & Draine, B. T. 2011, ApJ, 727, 35 
Dove, J. B. & Shull, J. M. 1994, ApJ, 430, 222 



10 



No extinction Constant extinction Prescott extinction 





1.0 
















0.8 


71% ; 






' 71% : 


71% ; 


Non-HII 
F(8um) 


ion 


0.6 








60% . 


.-• /' 62% . 


>10Myr 
F(Snm) 


Fract 


0.4 


> v 30% 


\ _ 




41% 


N - , ' ' 43% 


>100 Myr 
- F(8nm) 




0.2 


37% 






25% 




Diffuse 










F(Ha) 




0.0 








i ■ ■ i i ■ ■ ■ 


i ■ ■ ■ i ■ ■ ■ i ■ ■ ■ i ■ ■ ■ i ■ ■ ■ i ■ ■ ■ 




)/F(Ha) 


100 












Non-HII 

regions 

HII 

regions 














B 
3. 

GC 


10 
















1 















02468 10 12 02468 10 12 02468 10 12 
Radius (kpc) Radius (kpc) Radius (kpc) 



Figure 7. Differences between different extinction corrections. The left column shows the results when no extinction correction is applied. 
The middle and right columns both apply a radially declining extinction correc tion to the obs e rved n on-Hll region Ha emission, but the 
middle column applies a uniform extinction cor rection to the Hll re gions from Sanchez ct al. (2011]) while the right column applies the 
radially declining extinction correction found in Prescott ct al. ( 2007). The colors and line styles are identical to Fig. [3] and noted in the 
legend. 
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Figure 8. Testing the effects of different emission measure cutoffs. 
The colors denote AEM=1, 1.5, 2, 4 cm -6 , in rainbow order, from 
red to blue. The panels show from top to bottom: the non-H II Ha 
fraction, the non-H II 8^m dust fraction, and the 8^tm dust fraction 
heated by >100 Myr old stars. 
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APPENDIX 



A MODEL OF PAH EXCITATION BY YOUNG AND OLD STELLAR POPULATIONS 

Stellar spectral synthesis 

Stellar spectra c orresp onding to specified age ranges are created using the spectral synthesis code of 
iBruzual fc Chariot! ()2003l ) with Padova 1994 evolutionary tracks and a Chabrier IMF. We test metallicities of so- 
lar and two sub-solar values (z=0.02, 0.008 and 0.004) and a range of star formation histories. The first SFH is 
constant and the second is an exponential m odel with a timescale of 5 Gyr. The third is a linear model, based upon 
equation 9 of iFraternali fc Tomassett i (2012) and the 7 they determined specifically for NGC 628. T he final SFH is 
based upon the stellar population reconstruction from panchromatic photometry in IZou et aTl (|2011[ ) . This SFH is 
shown in their fig. 5 and peaks at 2 Gyr ago. We normalize all the SFHs to have a current SFR of 0.7 M Q yr _1 
(|Kennicutt et al.l 1201 If ) . This normalization results in total stellar masses (integrating over SFH and accounting for 
stellar mass loss) of 9.6, 10.1, 10.0 and 9.5 in lo g(Mf7)) for the const ant, exponential, linear and 2 Gyr-peaked SFHs. 
The measured log(M ) for NGC 628 is only 9.6 ([Skibba et al.H2011l) . Thus, the constant and 2 Gyr-peaked SFHs are 
accurately tuned to NGC 628, while the linear and exponential models have about three times too much stellar mass. 
Still, they may represent broadly the shape of NGC 628's SFH (perhaps a period of elevated SFR at the present) so 
we keep them for the analysis. 



Dust absorption: non-ionizing radiation 

In order to estimate the amount of PAH emission due to old stellar populations, we create a simple model based 
on separating the absorption of stellar energy due to PAHs from that of the total population of dust grains. We 
first assume that the old stellar populations are mixed with the dust while the youngest stellar populations follow the 
starburst attenuation law, which is broadly equivalent to a foreground (clumpy) screen. We change from the foreground 
screen (with optical depth r| cr ) to the mixed dust geometry (with optical depth r™ IX ) assuming a timescale for young 
clusters to emerge from their natal clouds of 10 Myr. The following equation describes the fraction of absorbed light 
as a function of the screen and mixed optical depths and the age of the stellar population (t, in Myr): 



fabs{T,t) 



- 1- 



l-e 



-t/w 



1 -e" 



(Al) 



The variation of optical depth with wavelength is computed by the ratio of the absorption cross-section of all dust 
(including PAHs), Ktot(A), normalized to that at V-band as follows: 



r(A) 



«tot(A) 
Ktot(V) 



Ty. 



(A2) 



We no te that this approach to the attenuation of an integrated stellar population is similar to that of Ida Cunha et aTl 
(2008), with an important difference being that we use a dust-model based cross-section instead of a simple power- law 
effective absorption curve. 

With Equations Al and A2 determining / a b s as a function of stellar population age and wavelength for given r v cr 
and Ty 11 ' 1 , we may then calculate the TIR luminosity for stellar populations with ages between t\ and ti as: 



^tir 



dt 



dX^-f a bs(t, A). 



(A3) 
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Figure 9. Left: The fraction of absorbed light versus wavelength for various age stellar populations, ranging from 1 Myr (blue) to 100 Myr 
(red). For these absorbed light fractions, we have assumed a T D m v =0.69 and a r youngi v =0.5, appropriate for the center of NGC 628. 
While these optical depths might seem backwards (higher optical depth towards older population), the different dust geometries assumed 
mean that the absorbed fraction is still higher towards the young population. 



For the PAH luminosity, we apply the ratio between the PAH-only cross-section to the total dust cross-section: 

£p AH = f 2 dt r^^^/ abs( i,A). 

Jti JXi dt «tot(A) 



Dust absorption: ionizing radiation 

For ionizing photons, hydrogen atoms present a major source of opacity in addition to dust. We assume that 
Hn regions are radiation bounded and that a fraction of ionizing photons, /ion, i onize hydrogen. The remainder are 
absor bed by dust within the H II region, which is assumed to be devoid of PAHs (e.g. lGiard et al.lll994UChurchwell et al.l 
2006). About 70% of ionizations produce a Lya photon upon recombination (valid for electron densities, n c < 10 3 
cm' 3 . lD"rainell20lTI ). which is then repeatedly scattered. After several scatterings, most of these Lya photons will be 

/esc, escape. Let Qo be the rate of ionizing photons in s" 1 and Ir be the 
erg). If the average energy of an ionizing photon is approximately 1.2/h, 



absorbed by dust, while a small fraction, 
ionization energy of hydrogen (2.179 x 10 
then the TIR luminosity from within the Hn region from ionizing photons is: 



L 



TIR, Hll(ionizing) 



1.2(l-/ ion )+0.7*-/ ion (l-/ csc ) 



(A5) 



Note the factor 3/4 gives the energy of the Lya transition by the Rydberg formula. 

H II regions are surrounded by dusty PDRs. The line emission from any escaping Lya as well as from optical hydrogen 
lines and metallic cooling lines should also be included here, as this is essentially reprocessed ionizing radiation. We 
assume that the PDR has a covering fraction (</>pdr) of 0.7. For the optical recombination emission, we know that the 
30% of recombinations that do not produce a Lya photon will produce a two-photon decay from the 2s state. Before 
decaying to the 2s state, the electron must have already lost 0.25/h through recombination lines. For the metallic 
cooling lines, we assume they produce 0.1/h per hydrogen ionization. Combining all of these terms, we have: 



LtIR, PDR(ionizing) — Qo^H/ion^PDR 



/esc + 0.3 * - + 0.2E 



o.i: 



PAHs do exist within the PDR, so we also have: 



LpAH, PDR(ionizing) — Qo^H/ion'/'PDR 



tpAH(Lya) 
retot(Lya) 



/esc 



«PAH(opt) 
Ktot(opt) 



0.575 



(A6) 



(A7) 



Because the dust within the Hn regions is assumed to be PAH- free, the ionizing radiation does not have a large 
contribution to the PAH emission. As shown in the left panel of Fig. [TU] (dotted black line), the escaping recombination 
and cooling radiation contributes less than 2% to the total PAH emission (assuming a constant SFR 10 Gyr-old 
population), even with a high /i on of 0.9. The TIR emission from intercepted ionizing photon light (either directly or 
through lines) is approximately 10% as seen in the right panel of Fig. [10} We note that thi s fraction agrees with the 
5-10% found for FIR emission from Hn regions within the Milky Way ([Sodroski et al.lll997[ ). 
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Figure 10. The fraction of the total PAH (left) and TIR (right) emission contributed by ionizing photons as a function of the fraction of 
ionizing photons that actually ionize hydrogen. Different colored lines represent the fraction for different age stellar populations. The solid 
black line represents the fraction for a continuous SFR population up to 100 Myr while the dotted line does the same up to 10 Gyr. 



Model results 

Assuming a value of f- lon = 0.7 and combining the non-ionizing and ionizing contributions we obtain both the TIR 
and PAH luminosity for young stars we have: 

/■100 Myr <.10,000A ,r 

Lrm, young = I.OIQo/h + / dt „ dX— ±f ahB (t, A) (A8) 

JO J912A at 

and 

/•100 Myr /.lO.OOOA it 

ipAH, young = 0.057Q /h + / dt . dX x PAH \ ( / abs (t, A) . (A9) 

JO J912A dt KTIR(A) 

The old population luminosities may be calculated using Equations A3 and A4 with the wavelength range 912-10000A 
and the stellar population age range 100 Myr- 10 Gyr. 

The choices for T v cr and t™ ix used in Eqn. Al should be set to match NGC 628. From Section 3.2, we have 
A(Ha)— 1.03 — 0.53.R/R25 for the H ii regions. Assuming the continuum is less attenuated than the Ha emission by 
the factor 0.44 (|Calzetti et alj|i"994l ) and transforming to V-band and to optical depth instead of magnitudes, we have 
Ty cr = 0.51 — O.26-R/R25. For the older populations, we again adopt the prescription for the disk regions from Section 3.2 
(linearly declining with radius), but this time use the full attenuation through the disk (and neglecting Milky Way 
foreground extinction). As the stellar photons exciting PAHs are traveling in all directions, the effective dust optical 
depth for absorption of these photons is likely t o be comparable to the full-thickness optical depth, hence adopting 
that value. Using the dust surface densities from lAniano et all (|2012f ) and converting to optical depth, we parametrize 
Ty lx = 1.20 — 0.92i?/R,25. (The different geometry assumed for young and old regions means that the higher r™ lx 
values give similar or lower absorbed fractions of stellar radiation than the lower Ty cr .) These parameterizations of 
t™ IX and Ty cr with radius result in Fig.fTTJ which shows the radial behavior of the fractions of PAH and TIR luminosity 
contributed by old stars (top left and right, respectively), the ratio of TIR to stellar luminosity (bottom left) and the 
ratio of PAH to TIR luminosity (bottom right). Fig. 11 also exhibits the effect of different star formation histories, 
with four different SFHs shown as different colored lines. 

Observationally, the TIR to stellar ratio is determined to be 0.33 in lSkibba et all (|2011[ ) for NGC 628. We plot this 
value as a horizontal dotted line in the lower left panel of Fig. [TT] The linear and exponential SFHs are close to this 
value at all radii, while the constant and Zou SFHs both predict too much infrared per unit stellar luminosity. All 
SFHs and choices have a L(PAH) to L(TIR) ratio of approximately 0.3. This is larger than the typical observed 
value, but the way that PAH emission is usually measured does not include the continuum emission under the PAH 
features nor any continuum emission emitted at longer (> 20/xm) wavelengths. Furthermore, we did not use a dust 
cross-section perfectly tuned to NGC 628, the observationally-determined fractional PAH content (<7pah) is actually 
about 80% lower than in the dust model we used, thus the PAH to TIR ratio should also decrease by this amount. 

As for the PAH luminosity from old stars, the linear and exponential SFHs are quite close to the global value we 
derive, although they do not follow the specific shape from our Fig. [3] We take this as an indication that it is indeed 
possible for approximately 30% of the PAH luminosity to be powered by old stars. The upper right hand panel of 
Fig. HU shows the corresponding plot for the TIR. While the cross-section of the bigger grains is weighted more towards 
longer wavelengths, the contribution from within the H 11 regions partly makes up for this lower sensitivity to young 
stars and the fraction of TIR emission from old stars is only about 10% more than for PAHs. 
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Figure 11. Upper left: The fraction of PAH luminosity contributed by old stars. Upper right: The fraction of TIR luminosity 
contributed by old stars. Lo wer left: The ratio of TIR to unattenuated stellar luminosity. The dotted line represents the global value for 
NGC 628 from Skibba et al. (2011). Lower right: The ratio of PAH to TIR luminosity. In all plots, the solid lines assume that the old 
population attenuation decli nes proportionally t o that of the H II regions while the dashed lines assume the old population's attenuation 
exponentially declines as per White ct al. (2000). The different colors reflect the different SFHs as shown in the legend. 
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